Splitting dynamics of giant vortices in dilute Bose-Einstein condensates 



Pekko Kuopanportti^Q and Mikko Mottonen 1, 2 

'Department of Applied Physics/COMP, Aalto University, P.O. Box 15100, FI-00076 AALTO, Finland 
2 Low Temperature Laboratory, Aalto University, P.O. Box 13500, FI-00076 AALTO, Finland 

(Dated: April 2, 2010) 

We study the splitting of multiply quantized vortices with large quantum numbers in dilute 
nonrotated Bose-Einstein condensates in the zero-temperature limit. The splitting is observed to 
result in vortex-free condensate fragments which are separated by vortex sheets. The number of 
these fragments is found to be equal to the angular-momentum quantum number of the Bogoliubov 
excitation mode responsible for the splitting, although the formulation of the fragments cannot be 
described by small-amplitude excitations. Thus, the realization of an isolated giant vortex and the 
observation of its splitting would provide a means to directly relate the experimental data to discrete 
theoretical quantities. 

PACS numbers: 03.75.Kk, 03.75.Lm, 67.85.De 



I. INTRODUCTION 

The creation and observation of quantized vortices in 
dilute Bose-Einstein condensates (BECs) of alkali-metal 
atoms in 1999 [l[ was an important demonstration of the 
supcrfluid properties of these systems. Since then, the 
study of vortices in BECs has nourished both theoreti- 
cally and experimentally 0. 

When a large amount of angular momentum is trans- 
ferred to a dilute BEC, it typically responds by nucleat- 
ing vortices. If a condensate in a harmonic trap is sub- 
jected to rapid rotation, large arrays of singly quantized 
vortices are observed [JQj- However, the condensate can 
also acquire high angular momentum through multiquan- 
tum vortices, for which the phase of the condensate order 
parameter winds an integer multiple k of 2ir. Multiquan- 
tum vortices have been created in dilute BECs by using 
a focused laser beam to remove atoms from the center 
of a rotating condensate 

by transferring angular 
momentum into the condensate by a Lagucrre-Gaussian 
laser beam [111 [l2j , and by a topological phase engineer- 
ing method |l3[ which utilizes the spin degree of freedom 
of the condensate [l3 and its coupling to an external 
magnetic field [H-EJ 

Several theoretical studies [22U27| have shown that 
states with a multiquantum vortex are typically dynam- 
ically unstable in harmonically trapped BECs. Dynami- 
cal instability is a peculiar feature of nonlinear dynamics, 
and in the case of multiquantum vortices, it causes the 
vortex to split into si ngle-qu antum vortices even in the 
absence of dissipation |28rl33| . On the other hand, vari- 
ous studies have also addressed different means to stabi- 
lize multiquantum vortices, e.g., by rotating the conden- 
sate in the presence of a plug potential 34| or an anhar- 
monic trapping potential [35l - l38| . Stable multiquantum 
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vortices have also been found to exist in two-component 

BECs (SUES. 

Recently, it has been suggested that vortices with arbi- 
trarily large winding numbers could be achieved in BECs 
with the so-called vortex pump, i.e., by cyclically pump- 
ing vorticity into the condensate [4l|. The technique is 
based on the topological phase engineering method: the 
spin degree of freedom of the condensate is controlled 
locally by adiabatically tuning external magnetic fields, 
which renders the system to acquire a fixed amount of 
vorticit y in each pumping cycle. In the cycle presented 
in Ref. [41[ , a homogeneous magnetic field is used along- 
side alternating quadrupole and hexapole magnetic fields. 
With this setup, the vortex pump can be operated both 
fully adiabatically and partly nonadiabatically. In the 
former case, an additional optical potential is needed to 
confine the condensate during the cycle. Furthermore, in 
order to stabilize the vortex and prevent it from split- 
ting prematurely, an optical plug potential can also be 
employed during the pumping stage. The plug can be 
realized by a focused laser beam as was done, e.g., in 
Ref. (Tfj| . The partly nonadiabatic operation, on the 
other hand, can be carried out without any optical po- 
tentials, but with the cost of lo sing a part of the atoms 
from the trap. Later, Xu et al. [42[ have investigated an 
alternative pumping cycle in which one of the multipolc 
magnetic fields is replaced by a second homogeneous field 
but which necessitates the use of an optical trap. 

In adiabatic vortex pumping and other methods based 
on topological phase imprinting, the condensate remains 
in its instantaneous eigenstate throughout the process, 
and the final state with a multiquantum vortex is very 
close to a stationary state. Therefore, the stability prop- 
erties of stationary multiquantum vortices become essen- 
tial for determining how large winding numbers can be 
reached with the pump. Dynamical instabilities and core 
sizes of giant vortices with large winding numbers were 
recently studied in Ref. |27j . The maximum strength 
of the dynamical instability of a K-quantum vortex was 
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observed to increase very slowly with n. On the other 
hand, the core size of the vortex, which partly determines 
the maximum adiabatic pumping speed, was found to in- 
crease roughly as y/R, for large k. The results of Ref. [27j 
suggest that giant vortices with extremely large winding 
numbers could be created by gradually speeding up the 
operation of the vortex pump. 

Motivated by the promising possibility of realizing iso- 
lated vortices with very large winding numbers, we in- 
vestigate the splitting dynamics and stabilization of such 
vortex states. We limit the scope of our study to pancake- 
shaped BECs in nonrotating harmonic traps and study 
how the giant vortices can be dynamically stabilized with 
the application of a plug potential in the center of the 
trap. In order to address the dynamics of vortex splitting, 
we numerically solve the temporal evolution of the sta- 
tionary vortex states subjected to small-amplitude per- 
turbations. Based on earlier studies for two- and four- 
quantum vortices [1^, |33| , we expect the splitting to be 
driven by dynamical instabilities and anticipate the split- 
ting patterns to reflect the symmetry properties of the 
excitation modes that trigger the splitting processes. 

The remainder of this article is organized as follows. 
In Sec. HH we present the zero-temperature mean-field 
theory of the condensate and describe how it is employed 
in the numerical calculations. Section IIIII presents our 
results on the stabilization and splitting of giant vortices. 
In Sec. IIV1 we discuss the main results of the work. 



\J h/mujz is the harmonic oscillator length in the axial 
direction. The order parameter is normalized such that 
||v]/|| 2 = J |\I>(r, cj>)\ 2 r drd0 = N, where N is the number 
of condensed atoms. The effective interaction strength 
<72D is related to the vacuum s-wave scattering length a 
by g2D = \f^y?ajma z and is assumed to be positive. 

Stationary states of the condensates satisfy the time- 
independent GP equation, which is obtained from Eq. ([T]) 
with the replacement ihdt — > fx, where fj, is the chemical 
potential. We choose the stationary states to be axisym- 
metric vortex states with a quantum number k, which 
implies that the order parameter can be written in the 
form 
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where f(r) is real and nonnegative. Small-amplitude 
oscillations about the stationary states can be studied 
in terms of the Bogoliubov quasiparticles J44|. The or- 
der parameter is decomposed in the form \T>(r,<fi,t) = 
exp(— [^(r, <j)) + x(r, <p, t)], where the small varia- 
tion x is written as 

X = J2 ["<zMe i(K+ ' ,)0 ~ ItV + v *(r)e^ K -^ + K* . (5) 



By substituting this decomposition into Eq. (TTJ) and lin- 
earizing with respect to x, we arrive at the Bogoliubov 
equations 



II. THEORY AND METHODS 

We restrict our studies to the zero-temperature limit 
and neglect possible finite-temperature effects. Ex- 
periments with dilute BECs can be routinely carried 
out at temperatures where this approximation is justi- 
fied [43| . We consider BECs consisting of bosonic atoms 
of mass m enclosed in nonrotating cylindrically symmet- 
ric harmonic traps and assume that the condensates are 
pancake-shaped, i.e., that the trapping frequencies in the 
axial and radial directions satisfy uj z 3> uj r . Under these 
assumptions, the order parameter \V of the condensate 
satisfies the Gross-Pitaevskii (GP) equation, 

»m*(r,0,t)= [H + .g 2D |*(r,<M)| 2 ]*(r,<M), (1) 

where the single-particle Hamiltonian is given by 

h 2 /„o 1„ 1 



U 



.d 2 r + -d r + ^dj)+V(r). 
2m \ r 



The potential V(r) is taken to have the form 
V(r) = \mu 2 .r 2 + V piug (r), 



(2) 



(3) 



where Vpi ug (r) denotes the possible optical plug poten- 
tial. In Eq. ([1]), the z dependence has been factored 
out as Vt'(r) = \l/(r, 4>) exp[— z 2 /2a 2 ]/ \Jira 2 , where a z = 



C K+lq 9f 2 {r) U%WU k ft! 1, (6) 



-gf\r) -C K ^ lq 



v q (r) 



v q {r) 



where 



Ck = -2^ 15 + r 9r ^) + V{r) [L + 29 f ' (7) 

The functions u q and v q are the quasiparticlc amplitudes 
corresponding to the eigcnfrcqucncy u q . Moreover, each 
excitation is characterized by an integer l q which deter- 
mines the angular momentum of the excitation mode 
with respect to the condensate. Without loss of gen- 
erality, we will assume that l q > (4f|. The Bogoliubov 
quasiparticle description is valid only if the oscillatory 
part x has a small norm compared with \V. 

The eigenfrequency spectrum {ui q } can be used to 
analyze the stability of a given stationary state. If 
the spectrum contains excitations with a positive norm 
/ [|it g | 2 — \vq\ 2 ] rdr but a negative eigenfrequency ui q , 
the stationary state is energetically unstable. Further- 
more, the stationary state is dynamically unstable if 
the quasiparticle spectrum contains at least one eigen- 
frequency with a positive imaginary part |46| . Occupa- 
tions of such complex-frequency modes initially increase 
exponentially in time [Eq. (|5j)], and thus small perturba- 
tions of a dynamically unstable stationary state typically 
lead to significant changes in its structure. The complex- 
frequency modes also quickly drive the system beyond 



3 



the linear regime of Eq. ([5]) , and hence the complete dy- 
namics must instead be described with Eq. ([1]). 

In the case of a multiquantum vortex, dynamical in- 
stability typically signifies that the vortex is unstable 
against splitting into singly quantized vortices. More- 
over, the splitting pattern of the vortex is related to the 
angular momentum of the complex-frequency mode that 
induces the splitting (33j . A mode with the angular- 
momentum quantum number l q corresponds to a defor- 
mation of the condensate density with Z g -fold rotational 
symmetry. Hence, if the splitting of a multiquantum vor- 
tex is driven by an excitation mode with the quantum 
number l q , the splitting should initially exhibit ^-fold 
symmetry. In general, complex frequencies have been 
found to exist only for l q > 2 [H, [H, H] . 

To address the dynamical stabilization of multiquan- 
tum vortices, we study the effect of applying a localized 
plug potential along the symmetry axis of the trap. For 
simplicity, we use a step potential to approximate a Gaus- 
sian laser beam with a sharp boundary, i.e., 



VpiugO) = AQ(R plus -r), 



(8) 



where 8 denotes the Heaviside step function. The depen- 
dence on the amplitude A is removed by assuming that 
the plug is sufficiently strong such that A 3> /i, in which 
case effectively all particles are removed from the volume 
occupied by the plug. In this limit of strong pinning, 
the difference between the steplike and Gaussian-shaped 
beams becomes negligible [34^ . 

Our aim is to investigate the dynamical stabilization 
and splitting patterns of giant vortices. To this end, we 
numerically study stationary K-quantum vortex states for 
k > 2 and their response to slight perturbations. First, 
we find the stationary states of the form of Eq. (j4j for 
different values of k, g^c-, and i? p i ug by solving the time- 
independent GP equation. Utilizing these states, we 
solve the respective Bogoliubov equations, Eq. ([6]), and 
search for excitation modes that correspond to cigenfrc- 
quencies with large imaginary parts. From the excitation 
spectra, we also infer the minimum plug widths i? p i ug 
required to dynamically stabilize the vortex states. Sec- 
ond, the obtained stationary states are perturbed slightly 
by adding either low-amplitude random noise or a small 
initial population of a complex-frequency quasiparticlc 
mode. The dynamics of the perturbed states are then 
solved by numerical integration of Eq. ((T|). In order to 
relate the dynamics to quasiparticle modes with different 
angular- momentum quantum numbers l q , we monitor the 
evolution of different angular-momentum eigenmodes by 
computing the projections 
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-^(r^W, (9) 



where k is an integer. Since the evolution described by 
Eq. (fTJ) is unitary, we have J2k 11^*11 2 = N at all times. 

In the numerical simulations, we use finite difference 
methods with square xy grids that have roughly 200 



points in each dimension. The stationary states are 
solved using relaxation methods, and Eq. (fT|) is numer- 
ically integrated with the Strang splitting scheme [47j . 
The discretized Bogoliubov equations are solved using 
the LAPACK numerical library implemented in MAT- 
LAB [48[. In order to have generally applicable results, 
we measure length in the units of the radial harmonic 
oscillator length a r = yjh/mu: r and time in units of 
l/uj r . By normalizing the dimensionlcss order param- 
eter to unity, the dimensionless interaction strength be- 
comes g = \/&KNa/a z . For example, if we use the values 
a = 4.7 nm and m = 1.44 x 10 _25 kg corresponding to 
87 Rb atoms [44| and set the radial and axial oscillator 
frequencies to u r — 2tt x 20 Hz and w z = 50 u> r [49[ , we 
obtain a r w 2.41 /zm and g ~ 0.069 A^. Therefore, the 
effective interaction strength for a 87 Rb condensate with 
TV ss 10 4 atoms is of the order of g ~ 700. 

We employ two different types of perturbations in the 
stationary K-quantum vortex states. First, we add a 
small initial population of a selected complex-frequency 
excitation mode q such that 



*(r,0,O) = f{r)e lK * + r) m u q (r)e 



-v*(r)e 



i(K—lg)tp 



(10) 



where < rj m -C 1 determines the initial mode 
population and the quasiparticle amplitudes are nor- 
malized such that 2 J \u q (r)\ 2 r dr = 2 J \v q (r)\ 2 r dr — 
J f 2 {r)rdr. Second, we study the effect of adding ran- 
dom noise in each lattice point (r^, <pj) of the discretized 
order parameter such that the perturbed state is given 
by 

VfafaO) = f(n)e^ +VnPi,jf(ri), (11) 

where < ?/ n <C 1 characterizes the strength of the noise 
and pij is a random variable chosen uniformly from the 
unit disk in the C plane. Different grid sizes are used to 
ensure that the effect of the noise is independent of dis- 
cretization. Since the random noise is expected to pop- 
ulate different excitation modes approximately equally, 
we anticipate the noise-induced splitting to be driven 
by the dominant complex- frequency mode, i.e., by the 
quasiparticlc mode q with the greatest imaginary part 
of the eigenfrequency ui q . Consequently, the splitting is 
expected to exhibit Zdom-fold rotational symmetry, where 
idom denotes the angular-momentum quantum number of 
the dominant complex-frequency mode. 



III. RESULTS 

We have solved the time-independent GP equation and 
Bogoliubov excitation spectra, Eq. ©J for different val- 
ues of the winding number k, dimensionlcss interaction 
strength g, and plug width i? p i ug . The parameters k and 
g were primarily chosen from the experimentally most 
relevant ranges 2 < k < 50 and < g < 1000, although 
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in some cases greater values of g have also been used. For 
each re, we have determined the smallest width i? p i ug of 
the plug potential, Eq. ([SJ, which is sufficient to dynam- 
ically stabilize the re-quantum vortex states in the whole 
interval < g < 1000. Furthermore, we have computed 
the temporal evolution of slightly perturbed stationary 
multiquantum vortices from Eq. (TTJ) for different values 
of re and g. The evolution was studied without the plug 
potential for two types of perturbations, i.e., the addition 
of either random noise or a particular complex- frequency 
eigenmodc. 

In the absence of a plug potential, vortex states with 
large enough winding numbers re were found to be possess 
dynamical instabilities at nearly all values of g: vortices 
with re > 9 proved to be dynamically unstable in the 
whole interval < g < 1000 [Fig. [Ha)]. On the other 
hand, dynamically stable vortex states with re > 9 were 
detected at larger interaction strengths [FigQJb)], but no 
regular behavior for these typically narrow stability re- 
gions could be inferred. Within the parameter ranges 
2 < re < 50 and < g < 1000, complex eigenfrequencies 
were found only for angular-momentum quantum num- 
bers l q satisfying 2 < l q < 13. For fixed values of re 
and g, there were typically several dynamical instabili- 
ties corresponding to different values of l q as shown in 
Fig. [IJb) . The number of different l q modes supporting 
complex eigenfrequencies tended to increase with re. 

The addition of a plug potential of the form of Eq. ([5]) 
was sufficient for dynamically stabilizing the giant vor- 
tex states. In Fig. [2ja), we plot, as a function of re, 
the limiting plug width i? p i ug at and above which the re- 
quantum vortex is dynamically stable in the whole inter- 
val < g < 1000. For simplicity, the results are presented 
only for the strong steplike beam with only one param- 
eter -Rpiug, but qualitatively similar behavior was also 
observed for the linewidth of a Gaussian-shaped beam 
with a fixed amplitude. The stabilizing value of i? p i ug 
was found to be significantly smaller than the radius of 
the vortex core, and thus the plug induced only minor 
changes in the condensate density profile [Fig. Hfb)] or 
energy. 

To study the different splitting patterns of giant vor- 
tices, we have numerically integrated Eq. ([1]) starting 
from a slightly perturbed stationary state with given val- 
ues of re and g. The optical plug was assumed to be ab- 
sent, Vpi U g = [Eq. ©]. These kind of initial conditions 
correspond, e.g., to the experimental situation in which 
the stabilizing plug is switched off adiabatically but much 
faster than it takes for the vortex to split. We also point 
out that the addition of a small-amplitude perturbation 
to a dynamically stable vortex state never resulted in 
splitting, implying that in this case, dynamical instabil- 
ity is a necessary condition for the splitting to occur in 
the absence of dissipation. 

When a particular complex-frequency mode was used 
as the initial perturbation, the simulations confirmed 
the expected correspondence between the splitting pat- 
tern and the symmetry of the driving mode: addition of 
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Figure 1: (Color online) (a) Maximum imaginary part of 
the eigenfrequencies, max, [|Im(u; 9 /a;r)|], in the (g, re) plane. 
Black areas correspond to vortex states that are dynamically 
stable, (b) Maximum imaginary part of the eigenfrequencies 
for the vortex state with k = 10 as a function of the inter- 
action strength g for fixed values of the angular-momentum 
quantum number l q given in the inset. For < g < 2000, the 
10-quantum vortex is dynamically unstable except for three 
narrow regions within the subinterval 1400 < g < 1950. The 
dimensionless interaction strength is given by g = \^8TrNa/a z , 
where N is the part icle num ber, a the s-wave scattering 
length, and a~ = ^/h/moj z the axial harmonic oscillator 
length. 



a complex-frequency excitation mode with the angular- 
momentum quantum number l q to a stationary multi- 
quantum vortex state resulted in the splitting of the vor- 
tex with / g -fold symmetry. 



Figures [3] and 0] illustrate 
the typical splitting patterns observed for l q = 2 and 
l q = 5, respectively. The populations of the most signif- 
icant angular- momentum eigenmodes [Eq. ©] arc 
also shown. The squared norms of the components Vk^ 
with k — re ± l q initially grow exponentially as predicted 
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Figure 2: (Color online) (a) Minimum plug width R p i us 
[Eq. |8])] required to stabilize vortices with a given winding 
number k in the whole interval < g < 1000. (b) Areal 
particle densities [/(r)] 2 of the stationary 10-quantum vortex 
states for g = 880 with (dashed curve) and without (solid 
curve) a stabilizing plug whose width i? p i ug = 1.50 a r is dis- 
played by the vertical line. 



by Eq. (J5]). The first appearance of single-quantum vor- 
tices in the density profile is associated with the onset of 
oscillatory behavior in the populations ||'P/- , I r || 2 , indicat- 
ing rapid mixing between a large number of modes with 
different angular momenta. 

When the dynamically unstable vortex states were per- 
turbed with random noise, numerical integration revealed 
a total of three different splitting patterns, referred to as 
linear, threefold, and fourfold splitting. The splitting 
patterns were independent of particular realizations of 
the noise such that only one type of splitting occurred for 
given values of k and g. The results are summarized in 
Fig. [5ja) , where the observed splitting patterns are com- 
pared with the values of the quantum number l q of the 
dominant complex-frequency modes. Figure [5jb) exem- 
plifies the three patterns for 20-quantum vortices. In lin- 
ear splitting, the multiquantum vortex disintegrates into 
a linear chain of single-quantum vortices. In threefold 
(fourfold) splitting, three (four) sheets of single-quantum 
vortices emerge symmetrically from the core of the gi- 
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Figure 3: (Color online) (a) Density of condensed atoms dur- 
ing the evolution of a 10-quantum vortex state after addition 
of a complex-frequency mode with quantum number l q = 2 
and amplitude ?7, n = 0.05. The interaction strength is set to 
g = 880 and the field of view in each panel is 16 a r x 16 a r . 
(b) Squared norms of the projections Vk^ for different values 
of k [Eq. ©] as functions of time for the evolution presented 
in (a). Solid, dashed, and dash-dotted curves correspond to 
k = 10, 12, and 8, respectively, which are the values yielding 
the largest norms. The dotted curve corresponds to the ex- 
ponential function exp[2Im(cj q )f] = exp(2 x 0.228 u) r t) given 
by the Bogoliubov equations. The time instants shown in (a) 
are indicated with vertical lines. The inset shows the same 
plot for t < 10/cj r on a semilogarithmic scale (the curve for 
k — 10 is not shown). 



ant vortex toward the surface of the cloud, leaving be- 
hind three (four) vortex-free condensate fragments which 
are separated from each other by the interstitial vortex 
sheets. Each of the fragments was observed to rotate 
about the trap center and about its own center. 

The linear, threefold, or fourfold splitting pattern ob- 
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Figure 4: (Color online) (a) Particle density of a 30-quantum 
vortex state with g — 200 after addition of a complex- 
frequency mode with l q — 5 and rj m = 0.05. The field of 
view is 18 a r x 18 a r . (b) Squared norms of the projections 
Vk^ for different values of k as functions of time for the 
evolution presented in (a). Solid, dashed, and dash-dotted 
curves correspond to k = 30, 35, and 25, respectively. The 
dotted curve corresponds to the exponential exp[2lm(uj q )t] = 
exp(2 x 0.236 uj r i) given by the Bogoliubov equations. The 
time instants shown in (a) are indicated with vertical lines. 
The inset shows the same plot for t < 10/o; r on a semiloga- 
rithmic scale (the curve for k = 30 is not shown). 



served after addition of random noise was related to the 
presence of a complex-frequency mode with the angular- 
momentum quantum number l q = 2, 3, or 4, respec- 
tively, in the excitation spectrum of the corresponding 
stationary state. Moreover, the noise-induced splitting- 
process was reproduced by using the respective complex- 
frequency excitation mode as the initial perturbation 
[Eq. pop]. For example, when the state with k = 10 and 
g = 880 was initially perturbed with random noise of 
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Figure 5: (Color online) Noise-induced splitting of giant vor- 
tices, (a) The observed splitting patterns and dominant 
complex- frequency modes of stationary multiquantum vor- 
tices in the (g, k) plane. The lines separate regions where 
^dom, the angular- momentum quantum number of the domi- 
nant complex-frequency mode, is constant and has the value 
indicated. The shaded area corresponds to alternating values 
of idom and regions of dynamical stability [cf. Fig.[]Ja)]. The 
bars, triangles, and squares indicate splitting patterns with 
2-, 3-, and 4-fold symmetries, respectively, as observed in the 
numerical integration of the GP equation after addition of 
random noise of strength g n = 0.05 [Eq. (|11|) ]. (b) Density 
profiles obtained for n = 20 and different values of g illus- 
trating the three types of splitting. The field of view in each 
panel is 18 a r x 18 a T . 



strength ?y n = 0.05, the evolution of the condensate den- 
sity was indistinguishable from the evolution depicted in 
Fig. El for the corresponding mode-induced splitting, ex- 
cept that it took until t = 31.0/u; r before the vortex core 
had become elongated [panel t — 11.0/w r in Fig. El^a)] . 
The longer waiting time of noise-induced splitting is ex- 
plained by the smaller initial population of the driving 
complex-frequency mode. Finally, we point out from 
Fig. [5ja) that the complex-frequency mode accounting 
for the observed splitting symmetry was not always the 
one with the largest imaginary part of oj q . 

We have also studied the effect of dissipation on the 
splitting dynamics by adding a phenomenological damp- 
ing term to Eq. ([!} through the substitution t — > 
(1 — i^i)t with < 7 <^ 1 and normaliz ing, the order 
parameter after each time step [2!| [3(| H2 ■ Figure [S] 
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Figure 6: (Color online) Density of condensed atoms during 
the evolution of a 30-quantum vortex state after addition of a 
complex-frequency mode with l q = 5 and under dissipation of 
strength 7 = 0.01 (see the discussion in the text). Except for 
the dissipation, the setup is identical to Fig. 2] At t = 0.0/o; r , 
the expectation value of the axial angular momentum was 
(L») = 30.0 Nh. 

depicts the evolution of a state identical to the one shown 
in Fig. Eta) but under dissipation of strength 7 = 0.01. 
Comparison of Figs.@|a) and[S]shows that dissipation has 
caused the condensate fragments to shift closer to each 
other, but otherwise the splitting process has remained 
unchanged. Similarly, the splitting patterns presented in 
Fig. [SJa) were found to remain valid with 7 = 0.01. 

IV. DISCUSSION 

In this paper, we have investigated the splitting dy- 
namics of multiply quantized vortices in harmonically 
trapped pancake-shaped BECs for different values of the 
winding number k and interaction strength g. The split- 
ting was studied both by directly computing the tempo- 
ral evolution of the condensate from the time-dependent 
GP equation and by employing the Bogoliubov equa- 
tions to predict the essential features of the splitting pro- 
cess. To trigger the splitting, the stationary vortices were 
perturbed by adding either random noise or a complex- 
frequency Bogoliubov excitation to the order parameter. 
In addition, we have studied the dynamical stabilization 
of the giant vortices by applying a steplike plug potential 
along the symmetry axis of the trap. 

Dynamical instabilities were found to govern the split- 
ting behavior. Here, their existence was a necessary con- 
dition for the splitting to occur, and vortices with large 
winding numbers turned out to be dynamically unstable 
at nearly all values of the interaction strength g. On the 
other hand, a given vortex state could be stabilized by a 
plug potential of a width somewhat smaller than the size 
of the vortex core. 

The splitting mechanisms were found to be insensitive 
to the exact form of the initial perturbations. Even the 
addition of random noise to a stationary state resulted in 
a splitting pattern that could be reproduced by perturb- 



ing the state with a single complex-frequency mode. Typ- 
ically, this excitation mode was the one with the maxi- 
mal imaginary part of the cigenfrequency, although some- 
times a complex-frequency mode with a smaller value of 
the angular-momentum quantum number l q was excited 
by the noise instead of the dominant mode [Fig. [SJa)] . 
In general, the addition of a complex-frequency mode q 
resulted in a splitting pattern with Z 9 -fold rotational sym- 
metry. Moreover, the splitting patterns were found to be 
robust against weak dissipation. 

The splitting of giant vortices with large winding num- 
bers proceeded through the creation of vortex sheets and 
resulted in separate, vortex-free domains of condensed 
atoms. The number of these domains was given by the 
angular-momentum quantum number l q of the excitation 
responsible for the splitting. The deformation of the gi- 
ant vortex core and the eventual fragmentation of the 
surrounding condensate region were accompanied with 
large oscillations in the populations of different angular- 
momentum eigenmodes, suggesting that the formation of 
vortex sheets is a strongly nonlinear phenomenon. 

Our results indicate that a giant vortex does not split 
by melting into a chaotic liquid of singly quantized vor- 
tices. Instead, the long-time dynamics preserves the rota- 
tional symmetry of the excitation mode which causes the 
splitting. Therefore, the vortex pump with the simplest 
pumping cycle may not be an efficient way to realize a 
strongly correlated vortex liquid phase [HJ as envisioned 
in Ref. [4l|. 

In the future, it would be interesting to study the split- 
ting in a fully three-dimensional geometry to find out 
if additional effects such as intertwining of the vortex 
cores [29M3H would lead to more chaotic dynamics. For 
example, the intertwining might cause the vortex-free do- 
mains to intertwine along the axial direction (53j . 

The experimental creation of a nearly pure giant vor- 
tex state in a pancake-shaped BEC, e.g., using the vor- 
tex pump, and the consequent observation of its splitting- 
would provide a test for the validity of the Bogoliubov 
quasiparticle description. The splitting of the conden- 
sate into distinct fragments is detectable with current 
imaging techniques, and would allow for the clear iden- 
tification of the angular-momentum quantum number of 
the driving excitation. In addition, it would be inter- 
esting to determine if the phase coherence between the 
separate condensate fragments could be lost due to the 
rapid phase variations caused by the interstitial vortices. 
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